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Abstract 

We investigate numerically the Multiple Quantum (MQ) NMR dynamics in systems of nuclear 
spins 1/2 coupled by the dipole-dipole interactions in the case of the pseudopure initial state. 
Simulations of the MQ NMR with the real molecular structures such as six dipolar-coupled proton 
spins of a benzene, hydroxyl proton chains in calcium hydroxyapatite and fluorine chains in calcium 
fluorapatite open the way to experimental NMR testing of the obtained results. It was found that 
multiple-spin correlations are created faster in such experiments than in the usual MQ NMR 
experiments and can be used for the investigation of many-spin dynamics of nuclear spins in solids 
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I. INTRODUCTION 

Method of multipe-quantum nuclear magnetic resonance (MQ NMR) in solids [lj] is a 
powerful tool for the investigations of structure and spin dynamics in solids Q]. Significant 
progress has been made in the simplification of ordinary NMR spectra {3]. Recently, higher- 
order MQ NMR experiments have been utilized to address one of the fundamental issues in 
quantum information processing, the preparation of the pseudopure state 4|. A quantum 
computation includes as an initial step preparation of the input state [5j. Conveniently, the 
quantum algorithms start with a pure ground state where populations of all states except 
the ground state are equal to zero. The realization of a pure state in a real quantum system, 
such as a spin system, requires extremely low temperatures and very high magnetic fields. 
To overcome this problem, a so-called "pseudopure" state was introduced p, uj ■ The density 
matrix of the spin system in this state can be partitioned into two parts. The first part of the 
matrix is a scaled unit matrix, and the second part corresponds to a pure state. The scaled 
unit matrix does not contribute to observables, and it is not changed by unitary evolution 
transformations. Therefore, the behavior of a system in the pseudopure state is exactly the 
same as it would be in the pure state. 

Recently, it has been suggested 8, 9} to consider the dipolar ordered state as the initial 



state for such experiments. As a result of using the dipolar ordered initial condition, many- 
spin clusters and correlations appear faster j^, 9] than in the ordinary MQ NMR experiments 
lj in solids and some peculiarities of multiple-quantum (MQ) dynamics can be investigated 
with these experiments. 

In the present article we consider MQ NMR dynamics when the initial condition is de- 
termined by the pseudopure state. Our motivation of performance of this work is defined 
first of all by that the many-spin correlations are created faster in such experiments then in 
standard MQ NMR methods and can be used for the investigation of many-spin dynamics 
of nuclear spins in solids. We consider dynamics of MQ NMR in systems prepared in the 
pseudopure state and calculate the resulting signal in MQ experiments. Note, that it is not 
necessary to make any changes in the scheme of the standard experiment in order to obtain 
non-zero signals of MQ coherences in the pseudopure state. Computer simulations of such 
experiments for six spin ring and eight spin linear chain are presented. 
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II. MQ NMR WITH THE INITIAL PSEUDOPURE STATE 



We consider a system of nuclear spins (s = 1/2) coupled by the dipole-dipole interaction 
(DDI) in a strong external magnetic field H which is directed along the axis z. The secular 



part of the DDI Hamiltonian lOj has the following form 

(1) 



H dz = Djk I jz h z ~ \ (It Ik + IJIt) 
j<k L 



where Dj^ = ^(1 — 3 cos 2 6^) is the coupling constant between spins j and k, 7 is the 

r jk 

gyromagnetic ratio, is the distance between spins j and k, 6jk is the angle between the 
internuclear vector rjl and the external magnetic field. Ij a is the projection of the angular 
spin momentum operator on the axis a (a = x, y, z); and I~ are the raising and lowering 
operators of spin j. 

There are many pulse sequences exciting MQ coherences during the preparation period. 
For a dipolar-coupled spin system, the multipulse sequence with an eight-pulse cycle [l| 
is known to be very efficient. The basic scheme of the standard MQ NMR experiments 
consists of four distinct periods of time: preparation, evolution, mixing and detection 
In the rotating reference frame [h]] the average Hamiltonian describing the MQ dynamics 
at the preparation period can be presented as |l| 

n MQ = n {2) + n { ~ 2 \ (2) 

where H^ 2 ^ = —\^j <k Djklflk- Then the evolution period without any pulses follows. 
The density matrix of the spin system, p(r), at the end of the preparation period is 

P (t)=U(t)p(0)U + (t), (3) 

where U(r) = exp(—ir(H^ 2 ' + H^ 2 ^)) is the evolution operator on the preparation period 
and p(0) is the initial density matrix of the system. The transfer of the information about 
MQ coherences to the observable magnetization occurs during the mixing period. Usually 
the thermodynamical equilibrium density matrix is used as the initial one for MQ NMR 
experiments. In the high temperature approximation the equilibrium density matrix takes 
a form: p(0) = I z . Here we consider MQ NMR dynamics with the initial pseudopure state 
when the density matrix can be described as: 

p(0) = |l) 1 ®|l) 2 ®....®|l) jV , (4) 



where |l) fc represents a k— th spin that is up and iV is the number of spins in the system. 
The method of creating the highly polarized spin states (4) in clusters of coupled spins 



was described previously [4j, . It is based on filtering multiple-quantum coherence of the 
highest order, followed by a time-reversal period and partial saturation. It is convenient 
to expand the density matrix, p(r), at the end of the preparation period of the MQ NMR 



experiment as 
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p(r) = J2Pn(r), (5) 
where the term PnfV) is responsible for the MQ coherence of the n-th order. One can find 

e- i5/ *Vn(T)e^ = e- mSt p n (r). (6) 

On the mixing period the spin system is irradiated with sequence of the pulses shifted on 
a 7r/2— phase regarding the pulses of the preparation period [lj]. As a result the average 
Hamiltonian describing evolution of the spin system on the mixing period changes a sign on 
opposite to the sign of the Hamiltonian (2), and the evolution operator, U(r), is replaced by 
the operator U + (t). Starting with the initial conditions (4), the density matrix, p(t), after 
the three periods of the standard MQ NMR experiment can be written as 

p(t) = U + (T)e~ iShh p(r)e i5tlh U(r), (7) 

where p(r) is the density matrix at the end of the preparation period according to Eq. 
(3) and t = 2r + ti, S is the frequency offset on the evolution period of the duration ti 
which is a result of applying the time proportional phase incrementation (TPPI) method 
[l]. The transfer of the information about MQ coherences to the observable magnetization 
occurs during the mixing period. The last unitary transformation in (7) with operator 
U(t) describes this period. The resulting signal after the mixing period, the longitudinal 
magnetization, M z (t) , is 

M z (t) = Tr{p(t)I z } (8) 

Using Eqs. (3) and (7) and the initial condition (4) it is convenient to present the formula 
(8) for the longitudinal magnetization, M z (t) , as follows 

M z (t) = Tr {e~ m ^p{r)e m ^p MQ {r)} (9) 

where 

P mq (t) = U(t)I z U+(t), (10) 



coincides with the density matrix at the end of the preparation period of the standard MQ 
NMR experiment with the thermodynamical equilibrium density matrix as the initial condi- 
tion at high temperature approximation The density matrix p MQ {r) can be represented 
in the following form 
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P 



MQ 



j) = Y,Pn Q ^) (11) 



where the term p* f< ^(r) is responsible for the MQ coherence of the n-th order and satisfies 
to the relationship of Eq. (6). By using Eqs. (6), (7) and (9) one can rewrite the expression 
for the observable signal in terms of the intensities of MQ coherences 

M z (t) = J2^' mSt Ur)- (12) 

n 

where J n {r) are the normalized intensities of MQ coherences are 

J n (r) = 2Tr{p n (r)p M n Q (r)}/N. (13) 
One can find from Eq. (13) that 

J n (r) = J- n (r). (14) 

It is well known that in the usual MQ NMR experiments, the sum of the intensities of all 
MQ coherences does not depend on time ^ n J n (r) = 1. This is also right for the MQ NMR 
in the pseudopure state. 



III. THE NUMERICAL ANALYSIS OF THE TIME EVOLUTION OF MQ CO- 
HERENCES IN THE PSEUDOPURE STATE 

We restrict ourselves to numerical simulations of MQ NMR dynamics of one-dimensional 
systems. For example, quasi-one- dimensional six dipolar-coupled proton spins of a benzene 
molecule [4| and hydroxyl proton chains in calcium hydroxyapatite Ca^{OH){P04)^ [3] 
and fluorine chains in calcium fluorapatite Ca^FlPO^s [13J] are suitable objects to study 
MQ dynamics in pseudopure states. The numerical calculations were performed for MQ 
NMR dynamics of linear chains and rings of 6 and 8 spins. The DDI coupling constant of 
the nearest neighbors is chosen to be Djt+i = D = Is -1 . Then the coupling constants of 



spins j and k are D 



sin(7r/AQ 1 



for ring and D/\j — k\ 3 for chain, respectively. In order to 



sin(7r(j— k)/N) 

compare the results of the numerical simulations with the analogous ones for the ordinary 
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MQ NMR dynamics we introduce the normalized intensities of MQ coherences for the n 
order, J° rd {r) 

Jn d (r) = Tr [fl* (r)p M n Q (r)} /Tr {P z } . (15) 

The dependences of the intensities of MQ coherences on the dimensionless time, t = Dt , 
of the preparation period in spin ring containing six spins are presented in Fig. 1 for the 
both initial states, equilibrium (a) and pseudopure (b). One can compare the intensities of 
MQ coherences in the suggested experiment with the standard ones. It is evident that the 
suggested method can be considered as a useful addition to the standard MQ NMR method. 
Fig. 1 demonstrates that the suggested method yields the intensities of MQ coherences of 
the fourth and six orders which several times higher than the analogous coherence in the 
standard MQ NMR experiment. It is clear from the insets of Fig. 1 that MQ coherence 
of the fourth and sixth orders, obtained with the initial condition (4), appear little earlier 
than in the usual MQ NMR in a ring of six spins. For example, or the pseudopure state at 
time Dt = 2 the MQ intensity of the fourth order coherence is 0.05 (red dash line) while 
for the usually used initial condition up to Dt = 3 the intensity J4 is closed to zero (black 
solid line). A tendency of the faster growth of MQ coherences of high orders takes place also 
for the linear chain containing eight spins (Fig. 2). This peculiarity is connected with the 
initial pseudopure state (4). As a result, many-spin clusters and correlations connected with 
MQ coherences appear faster than in the standard MQ NMR with the initial condition, 
p(0) = I z . The numerical calculations confirm also the results obtained in the previous 
section. In particular, the growth of MQ coherences occurs in accordance with Eq. (14). 

One can see from Fig. 2 that the observable intensities of the fourth and sixth orders 
in the linear chain containing eight spins can be negative in contrast to the ordinary MQ 
NMR experiments at high temperatures . At the same time, it was shown [15j 



that the intensities of MQ coherences can be negative in the standard NMR experiments 
at low temperatures and in the dipolar ordered state at high temperatures. In fact, in MQ 
NMR experiments the observable quantity is the longitudinal magnetization modulated by 
rf pulses. The distinct frequency components of the magnetization can have an arbitrary 
sign. 
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IV. CONCLUSIONS 



The MQ NMR method for the detection of MQ coherences starting from the pseudopure 
state is proposed. Investigations of MQ NMR dynamics in the pseudopure states can be 
considered as a supplementary method which complements the usual NMR in order to study 
structures and dynamical processes in solids. Many-spin clusters and many-spin correlations 
are created faster in such experiments than in the usual MQ NMR with the initial equilibrium 
condition without any correlation between the spins. In this paper we focused on simple 
one-dimensional examples but the physical picture obtained here is not limited to perform 
simulations and experiments in two-dimensional and three dimensional systems and open 
new ways for the study of many-spin systems. 
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Caption figures 

Fig. 1 (Color online) The time dependence of the intensities of the MQ coherences in 
a ring of six spins coupled by the DDI: (a) termal equilibrium initial state and (b) the 
pseudopure initial state (4). Black solid line: intensities of the zeroth order coherence J . 
Red dash line: intensities of the second order coherence J 2 . Green dot line: intensities of the 
fourth order coherence J 4 . Blue dash-dot line: intensities of the sixth order coherence J 6 . 
The insets show that the MQ coherences of fourth (a) and sixth (b) orders in the pseudopure 
state (red dash line) appear little earlier than in the usual MQ NMR (black solid line). 

Fig. 2 (Color online) The time dependence of the intensities of the MQ coherences in a 
linear chain of eight spins coupled by the DDI: (a) termal equilibrium initial state and (b) 
the pseudopure initial state (4). Black solid line : intensities of the zeroth order coherence Jo 
. Red dash line: intensities of the second order coherence J 2 • Green dot line: intensities of 
the fourth order coherence J 4 . Blue dash-dot line: intensities of the sixth order coherence 
J 6 . The insets show that the MQ coherences of fourth (a) and sixth (b) orders in the 
pseudopure state (red dash line) appear little earlier than in the usual MQ NMR (black 
solid line). 
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